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A mechanism of d x 2_ y 2-wa,ve superconductivity is proposed for the high-T c cuprates based on a 
spin texture with non-zero topological density induced by doped holes through Zhang-Rice singlet 
' formation. The pairing interaction arises from the magnetic Lorentz force like interaction between 

| the holes and the spin textures. The stability of the pairing state against the vortex-vortex inter- 

^"*^ . action and the Coulomb repulsion is examined. The mechanism suggests appearance of a p-wave 

pairing component by introducing anisotropy in the CuCh plane. 
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I. INTRODUCTION 



O 
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3-h ■ High-temperature superconductivity occurs in a moderately hole doped Mott insulatorA The symmetry of the 
Oh' Cooper pairs has been established to be d x 2_ y 2-w&vei 2 To explain the mechanism of d x 2_ y 2-w&ve superconductivity, 
several mechanisms have been proposed. Among others, the antiferromagnetic spin fluctuation theory suggest a 
d x 2_ y 2-wave pairing mechanism^ between d-orbital electrons at copper sites. Another d x 2 _ y 2 -wave pairing mechanism 
is proposed in the resonating valence bond (RVB) theory. 1 In the RVB theory, electrons are assumed to be spin- 
charge separated. In the slave-boson approach^ spinons, which carry the spin degrees of freedom, form a d x 2_ y 2-wa,ve 
T \ pairing state. This pairing state is associated with the short-range antiferromagnetic correlations fi- Combined with 
Bose-Einstein condensation of holons, which carry the electric charge, the pairing state of the spinons leads to a 
d x 2_ y 2-w&ve superconducting state. 

Both of these theories are based on a pairing state between copper site electrons which produce antiferromagnetic 
correlations. However, from the sign of the Hall coefficient^ the charge carriers seem to be doped holes. Optical 
conductivity measurements also support that the charge carriers are doped holes because the Drude weight in the 
optical conductivity! is proportional to the doped hole concentration. If we assume that the charge carriers are 
doped holes, then it is natural to expect that high-temperature superconductivity is based on a d x 2_ y 2-w&ve pairing 
• mechanism between doped holes. In this paper, we propose such a mechanism based on a picture for the doped holes. 
The picture is that each hole induces a spin texture that is characterized by a non- vanishing topological charge density 
Q\ . in the spin system through the suppression of the antiferromagnetic correlations. 

In Ref. 8 , a half-skyrmion spin textur o 9 : 10 formation is discussed in the single- hole-doped cuprate within the non- 
linear sigma model(NLSM) description 11 of the Heisenberg antiferromagnet. The formation of the half-skyrmion spin 
texture can be understood by using an analogy to a superfluid system: Let us consider a two-dimensional boson system 
which exhibits superfluidity in the ground state. The Bose-Einstein condensate is described by a Gross-Pitaevskii 
equation^ If we suppress the condensate at some point P, then we see that a vortex is formed around P from the 
analysis of the Gross-Pitaevskii equation. Note that the field described by the Gross-Pitaevskii equation is a classical 
. field. However, the non-condensed component near the vortex core is described by quantum bosons. In case of a 
£^ ' charged boson system, the condensate is locally suppressed by a magnetic flux, or a vortex. 

We can apply this result to the two-dimensional quantum Heisenberg antiferromagnet. In the Schwingcr boson mean 
field theory^ the Neel ordering is described by Bose-Einstein condensation of the Schwinger bosons^ If we assume 
that the doped hole forms a spin singlet with a localized spin at a copper site, then the condensate is suppressed at 
that site. Forming a Zhang-Rice singlet in the antiferromagnet plays a similar role of introducing a magnetic flux in 
a charged boson condensate. As in the case of a superfluid system, a vortex appears around the site. In terms of the 
spins, the vortex is a spin texture with the non-vanishing topological charge density^ that is given by 



<7( r ) = ^ n ( r ) ' l d M r ) x d y n{r)) , (1) 



where the unit vector n(r) describes the staggered component of the spins. Note that the spin at the core of the spin 
texture is a quantum spin analogous to the normal component near a vortex core in a superfluid, and so the spin 
can form a singlet state with the doped hole spin, which is the Zhang-Rice singlet. 16 Based on the Schwinger boson 
mean field theory, we find that the topological charge, which is obtained by integrating q(r) over a whole region, is 
quantized to ±l/2^i£ Because of this quantized value, the spin texture is called half-skyrmion. 
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As argued in Ref£, the dispersion of the half-skyrmion is in good agreement with the results of the angle-resolved 
photoemission spectroscopy (ARPES) in the undoped compounds^ The dispersion of the half-skyrmion excitation 
is the same as that of the quasiparticles in the 7r-flux phase with a dynamically induced mass42i211 In fact, the 
effective theory of the half-skyrmion is similar to that of the 7r-flux phase: the half-skyrmions are described by Dirac 
fermions with a U(l) gauge field interaction. Anomalously broad peaks observed in ARPES are associated with the 
self-trapping of the half-skyrmions due to the coupling to the longitudinal gauge fieldi^l 

In Ref£, we assume that there is Neel ordering. The presense of the Neel ordering is used for obtaining the 
quantized topological charge. In the absence of the Neel ordering, the topological charge is no more quantized because 
the antiferromagnetic correlations decay exponentially with a length scale given by the antiferromagnetic correlation 
length £a F . However, the system preserves the local relationship between the hole density and the topological charge 
density. We consider spin textures with such topological charge density and hereafter we call such a spin texture 
st aggered- spin- vor t ex . 

The rest of the paper is organized as follows: In Sec|H] we describe the model for the staggered-spin-vortex 
formation. The pairing interaction is discussed in Sec lIIII We solve the gap equation in Sec II VI and show that the 
stable pairing state has d x 2_ y 2-wave symmetry. Section is devoted to the conclusion. 



II. MODEL 



We assume that the holes and the localized spins are independent degrees of freedom. Such a model is, for instance, 
the spin-fcrmion model, which is derived from the d-p model by applying a canonical transformation. In such a model, 
the spin system is described by the antiferromagnetic Heisenberg model, 

IL ./^S. -S ; . (2) 

Here the summation is taken over the nearest neighbor sites and Si denote the spin moment at copper sites. As an 
effective theory of the spin system, we consider the CP 1 model j22 which is derived from the NLSM by representing 
n(x) in terms of the Schwinger bosons: n = , ~z a cr ff)Cr > z a : 



9 



S CP1 = - / d 3 xJ2 



4,1 



(3) 



where A sw is the antiferromagnetic spin wave excitation gap and c sw is the spin wave velocity. The coupling g is 
given by g = c sw / p s with p s = Z„J/A the spin stiffness constant. Z p is the renormalization factor due to quantum 
phase fluctuations. In the spin disordered phase, A sw ^ 0, and there is no antiferromagnetic long-range order. The 
antiferromagnetic correlations are characterized by the correlation length ^af = c sw /A sw . For the Schwinger bosons, 
the presence of the antiferromagnetic correlations is associated with the phase coherence among the Schwinger bosons. 
For the kinetic energy of the holes, we assume 

H h =^e k c\ a c ka , (4) 

ho 

where e/. = k 2 /2m with m the effective mass. This is an approximate form in the continuum. 

The coupling between the holes and the localized spins Sj has the form of the antiferromagnetic Kondo coupling: 

II K 4'X('/'J'«r (5) 



The parameter 3k is on the order of 1 ~ 3eV. Because of this strong Kondo coupling, holes suppress phase coherence 
of the Schwinger bosons locally through the Zhang-Rice singlet formation. The suppression of the phase coherence 
leads to a staggered-spin-vortex: 

Vxa = f^4( r )^( r ). (6) 

where V x a = d x a y — d y a x and V'+f—) i s the creation operator of the staggered-spin-vortex (anti-staggered-spin- 
vortex). The coefficient is determined by considering the limit of £af — * oo. In this limit, the spin texture is the 
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half-skyrmion. The gauge flux V x a is associated with the spin chiralityiSiSi Note that (V x a) = 0, because 
in the equilibrium there is an equal number of staggered spin vortices and anti-staggered spin vortices, and so 
(E s #|(r)^ s (r))=0. 

By integrating out the Schwinger bosons, we obtain the effective action of the gauge field a M . Since A sw ^ 0, the 
action takes the Maxwellian form. That is, the gauge field is massless. However, in the spin disordered phase, the 
antiferromagnetic correlation length £af is finite. Therefore, the gauge field propagator decays with a length scale of 
£af- To include this feature, we perform a duality mapping. We write z a {x) as 

z a {x) = py 2 exp(i6» + i9 v ). (7) 

Here the phase 6q is associated with the coherent motion of the Schwinger bosons and the phase 9 V is associated 
with the staggered-spin-vortices. Note that this form is applicable for the description outside the vortex cores. After 
performing a Stratonovich-Hubbard transformation in Eq.(PJ), we obtain the following Lagrangian density, 

c = j z Jl + u ll (d ll e + d IJ ,e v -a ll ). (8) 

We find d t j,J fl = Q by integrating out 9q. This constraint is satisfied by introducing a gauge field: 

1 

Wc define the vortex current by 



J ^ = t^c^ac^a- (9) 



3l = ^e„ vX d v d x 9 v . (10) 
In terms of the staggered-spin-vortex and anti-staggered-spin-vortex fields, the vortex density is given by 

p<-(r) = X>J(r)^(r). (11) 

S 

From the equation of the continuity, dtp + V • j v = 0, we find 

j" = E^ MV^-(Wl)V.], (12) 

where we have used the fact that the kinetic energy of the holes is described by Eq.QJ. The coupling between the 
gauge field and the vortex current is of the minimal coupling: 

C mt = -ifiAp. (13) 

Since Jo = V x A/2tt and Jo describes phase fluctuations of the Schwinger boson density, nonvanishing gauge flux 
V x A describes the suppression of the phase coherence. Contrary to the gauge flux V x a, the flux quantum for the 
Meissner phase of A is 2ir because staggered-spin-vortex carry a unit gauge charge. Therefore, the relation between 
the vortex density and the gauge flux is 

^](r)^ s (r) = -i-VxA. (14) 

s 

The gauge field is equivalent to the gauge field in the slave-boson mean field theory^ 



III. PAIRING INTERACTION 

The pairing interaction is obtained from Eas. (|14|) and l|13fl by eliminating the gauge field A^. We take the Coulomb 
gauge: 

V-A = 0. (15) 

In the momentum space, 

Aqx = — T^<?i Any = T-^q- (16) 
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From Ea. (|14f> we obtain 

A, = -27rJ>L^*+ 9l .. (17) 



The staggered-spin-vortex current operator is 

j 



q = Ei( k +§H<w (18) 



m 

k.s 



Substituting these equations into Ea. (|13f) . we obtain 

2ni sr-^ q x k' , + + 

hlt = ^ ~q^ S ^q^k's'^+q^ka, (19) 

q,s.s' 

where SI is the area of the system. Since we are interested in the Cooper pairing, we set k + k' + q = 0. Thus, we 
obtain 

2iri sr-^ k x k' , * t 

int = ~^ ^ Ik-k'P 5 ^k'sW-k',A-Ks'^ks. (20) 
k/k',s,s' 1 1 

By making replacement of s — > s' , s' — ► s, k' — ► — k', and k — ► — k, we obtain 

2ni y-^ kxk' * . 

hlt= ~^ ^ | k _ k /|2 a ^»^-fc',s'^-M'V'fc fl - (21) 



Therefore, we may write 



H ^ = -^h E rk^p (s + s ' ) ^'^- fc '^ v '- fe ' s ' v ' fes - (22) 

From this form, it is apparent that the interaction exists for either staggered-spin-vortex pairs or anti-staggered-spin- 
vortex pairs. The Hamiltonian is given by 



H = Y^H^ S \ (23) 



H (s) = E - ^ E ^rA^-.^-^ks. (24) 

k k#k' 1 1 

A similar pairing interaction was discussed in the composite fermion system at half-filled Landau levels In that 
system, the composite fermions are spinless fermion and there is no other index associated with internal symmetry. 
It was shownSSiSl that the pairing interaction leads to p x ± ip y -w&ve pairing state. The sign is determined by the 
direction of the applied magnetic field perpendicular to the system. This pairing state is consistent with numerical 
simulations 

Before going into the analysis of the gap equation derived from the Hamiltonian 12MI) and l|24fl. we show an intuitive 
view about the origin of the attractive interaction. According to Eq. Q14JI. the gauge flux V x A is produced by the 
doped holes. Suppose a hole passes with the Fermi velocity the region of the gauge flux created by another hole. 
Then, as schematically shown in Fig. ^ the motion of the former hole is equivalent to a charged particle motion under 
a magnetic field. Magnetic Lorentz force like interaction is induced between the two holes. Such an interaction leads 
to a chiral pairing state. In this pairing mechanism, the gap is scaled by the Fermi energy. Since holes carry either 
positive or negative gauge charge, we expect there are two chiral pairing states with opposite chiralities. The stable 
pairing state in the bulk turns out to be c? 2 ,2_ J; 2-wave pairing state as we shall show in the next section. 



FIG. 1: Pairing interaction between the doped holes. Shaded region represents the gauge flux V x A created by the hole sitting 
at the center. If another hole passes this region with the Fermi velocity as shown in the figure, a magnetic Lorentz force like 
interaction is induced between the two holes. 



IV. GAP EQUATION 



Now we apply the BCS mean field theory to Ea. (|23|l . First we consider staggered-spin- vortex pairs described by 
H^ + \ In the following, we set Xk = ipk+ and Xk 
ensemble and define the following mean field: 



V'fc+i to simplify notations. We consider the grand canonical 



where 



The mean field Hamiltonian is 



k'(^k) 



Ani k x k' 



kk' 



Ik — k'l 2 " 



(25) 



(26) 



(27) 



where = eft — ji with ji the chemical potential and the summation in the momentum space is taken over half of the 
Brillouin zone. The gap equation is given by 



.(+) - 



= -— y v kk , 

2rt ^ kk 



^tanh^i 



k'(#k) 



(28) 



with Ek 



IA 



(+)i 



A similar gap equation is analyzed in Ref^ in the context of the composite fermion pairing. We apply the analysis 
presented in Refp^ 5 . In order to solve the gap equation, we introduce the following AnsataSi: 



A { + } = A fc exp (-i#k) , 



where I is an integer. In Eq. 128|) the integral with respect to 9k> is reduced to 



d6- 



sm( 



cos 6 — r\ 



exp(i£6), 



(29) 



(30) 



with r] = (k 2 + k' 2 )/(2kk'). The function It(r}) is exactly calculated by setting z = exp(i6) and applying a contour 
integral in the complex plane. From Eq. (|28p. we see that if the interaction is attractive (repulsive) then the sign of 
h{r]) is positive (negative). For I > case, It(rf) > while I < case, h(rf) < 0, and Ii = o(t]) — 0. Therefore, the gap 
equation has solutions only for £ > 0. Thus, we obtain 



1 

2m 



dk' 



k'A k , 
Ey 



00 ,,,k'A k , (k 
dk 



Ey 



(31) 
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From the asymptotic forms in k — > oo and k — > 0, we take the following approximate form 



25 



_ J Ae F (k/k F ) e for fc < fcp , . 

k ~ { Ae F (k F /k) e for k > k F ' ld j 

where is the Fermi energy of the holes. The gap A is obtained from the following equation: 

dx ; + / dx ; = 1. (33) 

v/(x 2 - l) 2 + A 2 x 21 Ji v / (x 2 -l) 2 + A 2 x- 2e 

The largest gap A = 0.916 is obtained for the case of p-w&ve (£ = 1). The second largest gap A = 0.406 is obtained 
for the case of d-wave {£ = 2). The third largest gap is A = 0.264 for £ = 3. 

For anti-staggered-spin-vortex pairs, we may carry out the same analysis. Since the interaction term has opposite 
sign compared to the staggered-spin-vortex pair case, the gap function A^ ' has the following form: 

A<f 5 = A k exp (tfflk) , (34) 

with t > 0. 

From the above analyses, the linear combination of the gap functions is 

A k = (A k +) +A ( - ) )/2 = A k cos {ifo) . (35) 

The relative phase between A^ and A^ ' is set to be zero. (A non-zero phase can appear if there is a magnetic field.) 
For the case of £ = 1, the symmetry of the Cooper pair is p x . Such a state can be stablized only at the boundary of 
the sample in the square lattice. Meanwhile, £ = 2 case, that is, d x 2_ y 2-w&ve pairing state is stable in the bulk. Thus, 
we may conclude that the above pairing interaction leads to d x 2_ y 2-w&ve superconductivity. 

The above analysis is carried out for the staggered-spin-vortices. In terms of those fields, the hole spin states are 
implicit. For spin singlet pairing states, the spin states should be 

(xikXl), (36) 

and 

(xL fc xi>, (37) 

where a in \k denotes the hole spin state. The relative phase, which is arbitrary in the above analysis, is it for the 
spin-singlet pairing state: 

(xifeXl) = -(xl k xl)- (38) 

Now we consider the effect of the staggered-spin vortex-vortex interaction and the Coulomb interaction. Those 
repulsive interactions reduce the gap value. In the following we consider those effects separately. In order to evaluate 
the vortex-vortex interaction effect, we consider a two- vortex solution as follows: 

9 V = tan" 1 - + tan" 1 — — (39) 
x x — d 



The interaction energy between the vortices is 



J ry ^ f ,2 2 x r zx(x-d,y) 



V vv (d) = -Z- P j ^— . ( ^ )2+ - (40) 

Taking into account the fact that the range of the interaction is limited by the antiferromagnetic correlation length 
£af, we obtain 

V vv {d) -irJZpp In (41) 



In the momentum space, 



V vv (q) = ^^JZ p p, (42) 
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with the constraint, q > 2tt/^af- 

In the presence of the vortex- vortex interaction V vv (q), the gap equation is given by 



1 = 



dx — — 



x 



2£+l 



X" - 1) 

i-C 

dx 



x 2e 



A-21 



dx — == 

x 







V(x 2 - l) 2 + A 2 x 2t 1 - x 2 



f A 2 x~ 2i 
dx 



-21' 



i+C 



V(x 2 -l) 2 + A 2 /x 2£ a; 2 - 1 



where £ = 2tt /{Uf^Af) and the constant C v is 



(43) 
(44) 

We take C v ~ 0.5 as a reasonable value at x = 0.10 by setting pCa_f ~ 1> J/ e F ~ 2. For Z p , the value of Z p = 0.72 
is used which is estimated from quantum Monte Carlo simulations. 29 Figure [5] shows Uf^af dependence of the gap 
parameters for each £. The effect of the vortex-vortex interaction is large for strong antiferromagnetic correlations. 
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FIG. 2: The gap parameter A versus ^af^f /(2tt) 



The gap values are somewhat reduced by the vortex-vortex interaction. While the relevant parameter region would 
be ^af^f/(27t) = 1 ~ 2, around £AF&F/(27r) ~ 3, the f-wave gap becomes larger than the d-wave gap. Note that as 
an approximation, C v is fixed in the above calculation for simplicity. For more precise calculations, we require detail 
analysis of {.afIcf dependence of ~p and Z p . 

Another pairing mechanism based on merons, which presumably correspond to the spin texture with non-zero 
topological density, was discussed by Berciu and Johni22i2i From the analysis of an extended Hubbard model, where 
a nearest-neighbor Coulomb repulsion is included, it was suggested that meron-antimeron pairs lead to d-wave super- 
conductivity. In Ref^fi^i, the pairing interaction between merons and antimerons comes from the vortex-antivortex 
interaction, and the pairing interaction (11'2[) is not considered. Although the vortex-antivortex interaction is attractive 
interaction, the analysis of the gap equation for the vortex-antivortex interaction within our Ansatz shows that there 
is no pairing state for ^af^f I '(27r) < 13. In this analysis, the coefficient C v is fixed. More precise analysis requires 
£af dependence of C v . However, C v is expected to be an increasing function with respect to £,af, and the increase of 
C v leads to reduction of the gap values for large £af- 

Now we consider the effect of the Coulomb interaction. The Coulomb interaction between the holes is 



(45) 
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with e the dielectric constant. The gap equation with the Coulomb interaction is 



.21+1 



1 = dx 
Jo 

-c e 



x 2l+l/2 



-l-2i 



^{x 2 - l) 2 + A*x~ ./ 1 



dx 



dx 



Ji 



^{x 2 - l) 2 + A 2 x~ 2e 



^(x 2 - l) 2 + A 2 x 2t \ 2a; 



x 2 + l 



dx 



J( 



^(x 2 -l) 2 + A 2 /x 2e \ 2x 



x A + l 



where 



Jo VV - cos < 



a = 



Figure shows 1/e dependence of the gap A for each t. 
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FIG. 3: The gap A versus 1/e. 



V. CONCLUSION 



In this paper, we have proposed a mechanism of <i :E 2_ y 2-wave superconductivity based on a spin texture with 
non-zero topological charge density. The spin texture formation is based on the Zhang-Rice singlet formation in the 
background of the antiferromagnetic correlations. In terms of a gauge field that describes antiferromagnetic spin 
correlations, the spin texture is described by a gauge flux. The interaction between the flux and the hole current 
induces the magnetic Lorentz force like interaction between the holes. Such an interaction leads to chiral pairing 
states with opposite chiralities. It turns out that the stable pairing state in the bulk is the d r 2_ a 2-wave pairing state. 
The stability of the pairing state is examined against the vortex- vortex interaction and the Coulomb interaction. 

In our pairing mechanism, the p-wave state is unstable in the bulk. However, this state can be stabilized in the 
presence of anisotropy. Since the p-wave gap is much larger than the d-wave gap, we expect enhancement of the 
superconducting transition temperature if we can induce that component. If there is the p-wave component, then the 
parity is broken locally. While the time-reversal symmetry is not broken unless one chiral pairing state is suppressed. 
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